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Abstract 

The master differential equations in the external square momentum p 2 for the 
master integrals of the two-loop sunrise graph, in n-continuous dimensions and for 
arbitrary values of the internal masses, are derived. The equations are then used 
for working out the values at p 2 = and the expansions in p 2 at p 2 = 0, in (n — 4) 
at n — » 4 limit and in 1/p 2 for large values of p 2 . 



PACS ll.10.-z Field theory 

PACS ll.10.Kk Field theories in dimensions other than four 
PACS 11.15.Bt General properties of perturbation theory 
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1 Introduction. 



In a previous paper it was shown that a linear system of first order differential equa- 
tions in the external Mandelstam variables, to be referred to from now on as master 
equations, can in general be written for the master integrals of any Feynman graph. In 
this paper the master equations are written for the 2-loop sunrise self-mass graph ampli- 
tudes in n continuous dimensions and for arbitrary values of the masses as a function of 
the square momentum transfer p 2 . All momenta are Euclidean, so that p 2 > corresponds 
to spacelike p^. The master equations can be used for obtaining valuable information on 
the properties of the master integrals, for working out virtually any desired expansion, for 
obtaining closed analytic results (when possible) and, in any case, for the actual numerical 
calculation of the amplitudes. 

In the various sections of this paper we will show 

• how to use the integration by part identities , [|J for writing the master equations; 

• how to obtain the (analytic) values of the master integrals at p 2 = 0, in the n — > 4 
limit; 

• how to expand the master integrals in p 2 at p 2 = 0; 

• how to expand the master integrals in (n — 4) for arbitrary p 2 ; 

• how to solve analytically the master equations when two masses are equal to zero; 

• how to expand the master integrals in 1/p 2 for large values of p 2 and arbitrary values 
of the masses. 

The appendix contains a few ready-to-use recurrence formulae, obtained by solving 
the integration by parts identities for the considered graph ||. 

All the algebra needed to generate the identities, the equations and the considered 
expansions has been processed by means of the computer program FORM Q by J. Ver- 
maseren. 

While we deal with the 2 loop case only, it should be apparent from the discussion 
that some suitable extension of the results may apply also to more general cases. 
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2 The master integrals and their differential equa- 
tions. 



The 2-loop sunrise self-mass graph is shown in Fig.l. 




Figure 1: The 2-loop sunrise self- mass graph. 

The related amplitudes, in n continuous dimensions and fully Euclidean variables, can 
be written as 



A(n, m\, ml, ml, p 2 , -an, -a 2 , -a 3 , (3 X , (3 2 ) = 

d n k x r d n k 2 {p-k 1 ) fh (jp-k 2 ) fh 



(2tt)"- 2 7 (2tt)"- 2 {k\ +m\)^{kl + m|) Q 2((p- h - k 2 ) 2 + m|)°3 



(1) 



It is by now well known that, as a consequence of the integration by parts identities 
0, all the above amplitudes can be expressed in terms of four master integrals ||, which 
can be taken as 

Fo(n,m\,ml,ml,p 2 ) = A(n, m\, m\, m\,p 2 , -1, -1, -1, 0, 0) , 
F x {n,m\,m\,m\,p 2 ) = A(n, ml, m\, m\, p 2 , -2, -1, -1, 0, 0) , 
F 2 {n,m\,m\,m\,p 2 ) = A(n, m\, m\, m 2 3 ,p 2 , -1, -2, -1, 0, 0) , 
F 3 (n,ml,m 2 2 ,m 2 3 ,p 2 ) = A(n, m\, m\, m\,p 2 , -1, -1, -2, 0, 0) . (2) 

It follows from the definitions that the last three master integrals are the mass derivatives 
of the first, i.e. for * = 1,2,3 

Fi{n,m\,m\,m\,p 2 ) = _ ^~2 f o( n ' m i' m 2' m 3.J' 2 ) ■ ( 3 ) 
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The simplest way for deriving the required differential equations for self-mass amplitudes, 
which depend on a single Mandelstam variable only, consists in writing for them the 
familiar scaling equation (which can be obtained, for instance, by scaling in the definition 
Eq.(|T]) all the integration momenta by a same factor A and then differentiating with 
respect to A). 

For F (n, m 2 , m 2 , m 2 ,p 2 ) that procedure gives 







drrir 



2 d 
+ m3 « 



(n-3) ) F (n,ml,ml,ml,p 2 ) = , (4) 



( 2 ® + m 2 ^ f m s 
\ dp 2 1 dm 2 2 

i.e., on account of Eq.(|]), 

2 u 771 / 2 2 2 2\ / IP / 2 2 2 2\ , 2 tt> / 2 2 2 2\ 

p—F (n,m 1 ,m 2 ,m 3 ,p) = (n - 3)F (n,m 1 ,m 2 ,m 3 ,p ) + m 1 F 1 (n,m l ,m 2 ,m 3 ,p ) 

i 2 771 / 2 2 2 2 \ i 2p / 2 2 2 2\ /r-\ 

+ m 2 F 2 {n,m 1 ,m 2 ,m 3 ,p ) + m 3 F 3 (n, m x , m 2 , m 3 , p ) . (5) 
A similar equation holds for each of F^n, mf, m 2 , m 3 ,p 2 ), % = 1, 2, 3 

(n-4) Fi(n,m 2 ,m2,m3,p 2 ) = . (6) 



9p 2 1 9m 2 



, d 9 d 
+ m 2 



dm 2 



When carrying out the mass derivatives, new integrals appear, like for instance those 
corresponding to «i = 3, a 2 = a 3 = 1 and (3% = /3 2 = in the notation of Eq. ([[]). All 
those integrals, according to the formulas of the Appendix, can be expressed in terms of 
the master integrals. The result can be written as 



^ u 771 / 2 2 22 

V — F i{ n , rn 1 ,m 2 , m 3 ,p 



n 



dp 2 



fc 771 / 2 2 2 2 



+ 



D(m 2 , m 2 , m 3 ,p 2 



+ 



D(m 2 , m^, rn 2 ,p 2 



r-i / 2 2 2 2\ ^ 7-t / 2 2 2 2 

P ifi {m 1 ,m 2 , m 3 , p ) — - — F (n, m l5 m 2 , m 3 ,p 



, -n f 2 2 2 2\ T7 f 2 2 2 2\ 

+ P ij i(m 1 ,m 2 ,m 3 ,p )F 1 {n,m 1 ,m 2 ,m 3 ,p ) 

,0/2 2 2 2\ tt> / 2 2 2 2\ 

+ Pi,2{mi,m 2 , m 3 ,p )F 2 {n,m 1 ,m 2 ,m 3 ,p ) 

,0/2 2 2 2\t->/ 2 2 2 2\ 

+P ij3 (m 1 ,m 2 ,m 3 ,p )F 3 (n,m 1 ,m 2 ,m 3 ,p ) 
Qi 5 i(m 2 , m 2 , m 2 , p 2 )T(n } m 2 )T(n, m 2 ) 

+ Qi,2(ml, rn\, ml, p 2 )T(n, m 2 )T(n, m 2 ) 
+Qi,3{ m i > m 2> ml,p 2 )T(n, ml)T(n, ml 



(7) 



3 



where D(m\, m|, m 2 ,p 2 ), the Pij(m 2 , m 2 ., mj,p 2 ) and the Qij(m\ ) m\ ) m\ ) p 2 ) are poly- 
nomials in mf,m 2 , 77t| and p 2 , to be discussed in a moment, while T(n, m 2 ) is defined 
by 



T(n, m 2 



1 



n-2 



(2vr) n - 2 P + m 2 
C(n) being the following function of n 

The expansion of Eq. @ in n around n 

C(n) = 1 - (n - 4) 
.1 



(n-2)(n-4) 



C{n) , 



r 3 



(8) 



(9) 



4 is 



log(2V^r) - -7e 



+ (^-4) 2 - 
+ G((n-4) : 



log 2 (2v^) + i[^ + 7 l 



- 7£;log(2\/7r) 



(10) 



where 7# is the Euler-Mascheroni constant. The expansion, which is given only to make 
easier the comparisons with other results in the literature, is however never needed in 
practice, despite the frequent presence of l/(n — 4) factors, as in Eq.(|S|). Indeed, there 
is always an overall factor C{n) for each loop; in the calculation of any physical (and 
therefore finite) quantity, the singularities, which are of the form C 2 (n)/(n — 4) 2 etc., 
cancel out, and in the remaining finite part it is sufficient to use the value at n — 4, which 
is simply 

C(4) = 1 . 
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The explicit expression of the polynomial D{m\, m 2 ., m 2 ,p 2 ) is 
D(m 2 ,m 2 ,m 3 ,p 2 ) = (p 2 f + 4(m 2 + m 2 . + ml)(p 2 ) 3 

+ 2 (?>m\ + 3m 2 + 3m 3 + 2m 2 m 2 + 2m 2 m 3 + 2m 2 m 3 ) (p 
+ A{m\ + m\ + f»3 — m\m\ — m\m\ — m\m\ — m 2 m 

+ 10m 2 m 2 m 3 )p 2 + i? 4 (m 2 , m 2 , m|) 
= (p 2 + (mi + m 2 + m 3 ) 2 ) (p 2 + (m x + m 2 - m 3 ) 2 ) 

■ (p 2 + (mi - m 2 + m 3 ) 2 ) (p 2 + (mi - m 2 - m 3 ) 2 ) 



2\2 



m 2 m 3 



:i2) 



where 



d2/ 2 2 2\ 

R (mi,m 2 ,m 3 J 



4 i 4 i 4 22 22 o 2 2 

m 1 + m 2 + m 3 — 2m 1 m 2 — 2m 1 m 3 — 2m 2 m 3 
(mi + m 2 + m 3 )(mi + m 2 - m 3 ) 
•(mi — m 2 + m 3 )(mi — m 2 — m 3 ) . 



(13) 
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The polynomial Pij, for i — 1 are given by 

P lfi {ml,ml,ml,p 2 ) = -(p 2 ) 3 - (3m 2 + m 2 . + ml)(p 2 ) 2 

— (3m^ — m 2 — mj — 1m\m\ — 1m\m\ + lQm\m\)p 2 

/ 2 2 2\d2/2 2 2\ 

— (m 1 — m 2 — m 3 )K [m 1 , m 2 , m 3 ) , 

o / 2 2 2 2\ o/2 2 2 2\ 

P 2 , {m 1 ,m 2 ,m 3 ,p ) = P lfi {m 2 ,m 1: m 3 ,p ) , 

n / 2 2 2 2\ n / 2 2 2 2\ 

P 3:0 {m 1: m 2: m 3: p ) = P lfi (m 3 ,m 2 ,m 1: p ) , 



Px^m 2 ,™ 2 ,™ 2 ^ 2 ) = -(2m 2 + m 2 + m 2 )(p 2 ) 3 

— (6m^ + 3m 2 + 3mg + 3m 2 m 2 + ?>m\m\ + 2m 2 m 2 )(£> 2 ) 2 

— (6m^ + 3m 2 + 3mg — hm\m\ — 5m^m 2 — 4m 2 m 2 — Am\m\ 
— 3m\m\ — 3m\m\ + 40m 2 m 2 m 2 )j9 2 

— (2m^ + wij + mg — 3m 2 m 2 — ?>m\m\ — 2m\m\)R 2 (m\, m 2 , m^) , 

o / 2 2 2 2\ o/2 2 2 2\ 

P2,2{m 1 ,m 2 ,m 3 ,p ) = Pi,i{m 2 , m ± , m 3 ,p ) , 

o / 2 2 2 2\ o/2 2 2 2\ 

P3,3{m 1: m 2: m 3: p ) = P 1A (m 3 ,m 2 ,m 1: p ) , 



o / 2 2 2 2\ 

Pi,2{m 1: m 2: m 3: p ) 



= m 9 



3(P 2 ) 3 - (7m 2 + 5m 2 - 3m 2 ) (p 2 ) 



„2\2 



— (bm\ + wij — 7mg — 6m 2 m 2 + 2m 2 m 2 + lAm\m\)p 2 

I 2 2 2\ d2/ 2 2 2\ 

— (m x — m 2 — m 3 )K {m 1 ,m 2 , m 3 ) 



n I 2 2 2 2> 

P li3 (m 1 ,m 2 ,m 3 ,p ) 

n 1 2 2 2 2\ 

P 2< i(m 1 ,m 2 ,m 3 ,p ) 

n / 2 2 2 2\ 

P2,3ijn 1 ,m 2 ,m 3 ,p ) 

n I 2 2 2 2^ 

P^x{m Xl m 2l m 3l p ) 



r> I 2 2 2 2\ 

P 1:2 (m 1 ,m 3 ,m 2 ,p ) , 

o / 2 2 2 2\ 

Pi^m^m^m^p ) , 

n / 2 2 2 2\ 

P h2 (m 2 ,m 3 ,m 1: p ) , 

o / 2 2 2 2\ 

P li2 (m 3 ,m 1: m 2 ,p ) , 

o / 2 2 2 2\ 

Pi,2{m 3 ,m 2 ,m 1: p ) ; 



n / 2 2 2 2\ 

P 3j2 (m 1: m 2 ,m 3 ,p ) 
similarly, the Qij are given by 

Q hl (m 2 l7 ml,m 2 3 ,p 2 ) = Q 2)2 {m\,m 2 2 ,m\,p 2 ) = Q 3 , 3 (m 2 ,m 2 2 ,m 2 3 ,p 2 ) 

„2\2 



3(p 2 ) 2 + 2{m{ + m 2 2 + m 2 3 )p 2 - R 2 (m(, m 2 , m 3 ) 



Q 



I 2 2 2 2\ 

i,2K, m 2 ,m 3 ,p ) 



4mf 



2\2 



(p 2 ) 3 + (m 2 + m 2 + 3m 2 ) (p 2 ) 
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— (m 1 + m\ — 3m 3 ~ \Qm\m\ + 2m\m\ + 2m\m\)p 2 

-(m\ + m 2 2 -m 2 3 )R 2 (m 2 1} m 2 2} m 2 3 )\ , (19) 

zO / 2 2 2 2\ ,0 / 2 2 2 2\ 

Qi i3 (m 1 ,m 2 ,m 3 ,p ) = Q lj2 {m x ,m 3 ,m 2 ,p ) , 
Q 2 ,i(ml,m 2 2 ,ml } p 2 ) = Q 1>2 {m 2 2 , ml, ml,p 2 ) , 
Q 2)3 {m\,m\,m 2 3 ,p 2 ) = Q^ 2 {m 2 2 ,m 2 3 ,m\,p 2 ) , 

-O / 2 2 2 2\ ,0 / 2 2 2 2\ 

Q 3 ,i( m i> m 2> m 3^ ) = Qx,A m ^ m \i m frV ) , 

Q3,2( m i> m 2> m hp 2 ) = Ql,2( m h m 2> m i>P 2 ) ■ (2°) 

The values corresponding to the other indices can be obtained by a suitable permutation 
of the masses. 

The above polynomials satisfy the relations 

22,0/2 2 2 2\ 1/2 2 2\r>/2 2 2 2\ 

m 2 m 3 Q 1A (m 1 ,m 2 ,m 3 ,p ) = --r{m x - m 2 - m 3 )P lfl [m 1 ,m 2 ,m 3 ,p ) 

1 ( n / 2 2 2 2\ o/2 2 2 2\ r> I 2 2 2 2 

+ g ^i,i( m D m 2, m 3,P ) - Piflijn^m^m^p ) - Pi^m^m^m^p 

22,o/2 2 2 2\ 1/2 2, 2\ n / 2 2 2 2\ 

l/o/2 2 2 2\ o/2 2 2 2\ , o / 2 2 22 

- - (Pi,i(m 1; m 2 ,m 3 ,p ) - P h2 (m x ,m 2 ,m 3 ,p ) + Pi i3 {m x ,m 2 ,m 3 ,p 

22,0/2 2 2 2\ 1/2, 2 2\n / 2 2 2 2\ 

m x m 2 Q 1>3 {m x ,m 2 ,m 3 ,p ) = -(m x + m 2 - m 3 )Pi i0 (m 1 , m 2 , m 3 , p ) 

- - (Pi^m^m^m^p 2 ) + P^ 2 {m\,m 2 2 ,m\,p 2 ) - P li3 (ml,m 2 2 ,m 2 3 ,p 2 )^ , (21) 

and 

7-1/ 2 2 2 2\ I 2 2 2 2\ o / 2 2 2 2\ 

D(m x ,m 2 ,m 3 ,p ) = -(p - m x - m 2 - m 3 )P lfi {m x ,m 2 ,m 3 ,p ) 

T> I 2 2 2 2\ t> / 2 2 2 2\ ry / 2 2 2 2\ /oo\ 

- P hX {m x ,m 2 ,m 3 ,p ) - P Xt2 (m x , m 2 , m 3 , p ) - P lj3 {m x ,m 2 ,m 3 ,p ) . (22) 

The equations Eq.s(^,^) are the master differential equations referred to in the title of 
this paper. They are long, but their structure is simple: they express the derivatives with 
respect to p 2 of the four master integrals as a linear combination of the master integrals 
themselves and of known non-homogenous terms, the coefficients being ratios of known 
polynomials in the masses and p 2 . The known terms are just products of the 1-loop, 
1-denominator vacuum amplitude T(n,m 2 ) of Eq.(|8|), the simplest integral in the game 
which fixes the overall normalization. 
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All the coefficients of the linear combinations are equal to a simple factor, depending 
on the dimension n times a ratio of polynomials in p 2 and the masses. There is a unique 
denominator, the polynomial D(m\ 1 m\,m\ 1 p 2 ) of Eq. (|12j) ; note that its zeros in p 2 cor- 
respond to the threshold and the three pseudo-thresholds of the considered 2-loop sunrise 
graph. 

It should be apparent, from the derivation, that master equations of similar structure, 
(i.e. equations whose coefficients are ratios of polynomials) can be established for the 
master integrals of virtually any Feynman graph, for any number of loops; the size of the 
polynomials is expected to grow quickly with the number of the denominators and of the 
loops, but they are in any case finite polynomials, which can be conveniently dealt with by 
means of an algebraic program. Quite in general, the non-homogeneous terms correspond 
to amplitudes of related Feynman graphs in which one of the propagators is missing; from 
this point of view, the Z-loop sunrise amplitudes are the simplest /-loop amplitudes, in 
the sense that in their equations the non-homogeneous terms are just products of I — 1 
1-denominator vacuum amplitudes, Eq.(g). Sunrise amplitudes are in turn present as 
non-homogeneous terms in the equations for the amplitudes of graphs having one more 
propagator and so on. 

As it will be shown in the rest of the paper, the master equations contain virtually all 
the information needed for studying the analytic properties of the master integrals and for 
obtaining any expansion of interest; furthermore, after some minor work (the expansion in 
(n— 4) at n = 4), the differential equations can be used for the actual numerical evaluation 
of the master integrals, without explicit reference to the original integral representation 
Eq.(|T]), once the values of the integrals are known at some initial value of p 2 . Both the 
expansion in (n — 4) and the evaluation of the initial values (at p 2 = 0), can be carried out 
by exploiting the information contained in the master equations (which hold identically in 
n and are regular at p 2 = despite the appearance of an apparent kinematical singularity 
at that value of p 2 ). 

3 The value at p 2 = 0. 

From the definitions Eq.s(|l|,2), it is obvious that all the master integrals as well as all their 
derivatives are regular at p 2 = (which in the considered Euclidean region corresponds 
to p^ = 0). If we introduce the vacuum graph of Fig.2 
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Figure 2: The 2-loop sunrise vacuum graph, 
and the corresponding amplitude 

rrf 2 2 2\ / ^"^1 f d H k 2 1 , . 

V[n, mi ,m 2 ,m 3 )-J {2n)n _ 2 J {2n)n _ 2 ^ + w^(^+m^((,k 1 + k 2 ) 2 + m^) ' 1 ' 

the p 2 = values of the master integrals defined in Eq.s(§) can then be written as 

F o (n,ml,ml,ml,0) = V(n,ml,n%,Tn%) , (24) 

d 

Fi(n,mi,m2,m|,0) = —^—^V(n,ml,ml,m 2 3 ) . (25) 

Let us now consider the differential equations Eq.([5|) and Eq.([7p for, say, % = 3 at p 2 = 0. 
Due to the regularity of the derivatives, the left hand sides of both equations vanish; 
therefore, the right hand sides must also vanish at p 2 = 0. It turns out that both equations 
involve the same combination of Fi(n, m\ , m 2 , m\, 0) and Fzin, m 2 , m 2 , m 2 ,, 0), namely 
m 2 Fi(n, mf, m 2 , m 2 , 0) + m 2 F 2 (n, m 2 , m 2 , m|, 0); by eliminating that combination, one is 
left with a single equation involving only Fo(n, m\, m 2 ,, rn 2 , 0) and i^n, mf, m|, mf, 0), 
which is its derivative with respect to m^, according to Eq.fl2"5|) ; in terms of the above 
vacuum amplitude the equation reads 



d2/ 2 2 2\ u \rl 2 2 2\ 

R (m 1 ,m 2 ,m 3 )-^^V{n,m 1 ,m 2 ,m 3 ) 



(n — 3)(m 3 — m 1 — m 2 )V(n, m Xl m 2} m 3 ) 

, / n\ m 3 ~ m l m 2rr/ 2\mt 2\ 

+ {n- 2) — | T(n, m 1 )T(n, m 3 ) 



+ (n - 2) mi T(n, m 2 2 )T(n, m») 
- (n - 2)T(n, m\)T(n, m 2 2 ) , (26) 



where R 2 (m 2 ,m 2 ,m 3 ) is the same as in Eq. (|T3"D . Eq.fl2"6D is the master equation for the 
2-loop vacuum amplitude of Fig. (2) - this time in one of the masses (the differential 
equations in the other masses are immediately obtained by a suitable permutation of the 
masses themselves). It is also easy to see that the same equations can be obtained by 
writing down directly the integration by part identities for V(n, mf, m 2 , m 3 ), Eq.(^) and 
the related amplitudes. 
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Eg . (|26|) can be used for evaluating V(n, m\, m\, mf) once its "initial" value corre- 
sponding to a particular value of ml is given. Note that the equation itself can provide 
such a value. Indeed, the vacuum amplitude and all its derivatives are regular for non- 
vanishing values of the masses; by considering the particular value 777.3 = mi + m2, the 
factor R 2 (ml, m\, m|) on the l.h.s. of Eq. (EH) vanish, so that the r.h.s. must also vanish; 
that amounts to the relation 



V(n, m 1; 777,2, ( m i + m 2)' 



77 



2(ti - 3) 
/ T(n, ml 



T(n, ml) T(n, 77t|) 



mi 



m 2 



T(n, ml) \ T(n, (mi + m 2 ) 2 ) 



mi 



m 2 



(mi + m 2 ) 



• (27) 



Eq. (E7|) can be taken as the "initial condition" which determines completely the solution 
of Eq . (|26|) . Therefore, the differential equation Eq . (p6|) , a first order equation in a single 
unknown function, can be solved in quadrature using the initial value given by Eq.( 



For the explicit calculations it is convenient to expand the amplitude V(n, m\, m 



2' 



m 



around n = 4, where it is expected to develop polar singularities in (n — 4). More exactly, 
Eq.([8|) shows that T(n,m 2 ) has a single pole (n — 4), so that the non homogenous part 
of Eq. (p6|) and the initial value V(n, m\,m\, (mi + m 2 ) 2 ), Eq. (P7|) , both develop double 
poles; it follows that also V(n, ml, m\, m|) can have at most double poles. It can therefore 
be expanded as 



17/ 2 2 2\ 

V{n,mi,m 2 ,m 3 ) 



C 2 (n 




4) : 



-V^ 2 \ml,m\,m 



(n -4) 



+^°)(m 2 ,m 2 ,m 2 ) + C(77-4) 



\/(^(m 2 ,m 2 ,m 2 ) 



(2f 



where C(n) is given by Eq.(]|) and C 2 (n) has been treated as an overall factor according 
to the remarks preceding Eq . (|TTD . Correspondingly, Eq . fl26|) splits into the following set 
of equations for the residua and the finite part at n = 4 



R 2 (ml, ml, ml) 







dm\ 



2 \m\,m\,m\) = (m\ — ml — m\)V^ 2 \ml,m\,m 



(ml — ml) + m\(m\ 



m q 



(29) 



R 2 (m 2 i,m 2 2 ,ml) 







dm': 



■V^ 1 \ml,m\,m\) = (m\ — ml — m\)V^ ^(m^m^m 



+ (ml — ml — m\)V^ 2 '(m\,ml,m 
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, 1 f 2/ 2 2\ , 2/2 2\ 

+-{m 1 (m 1 -m 3 ) + m 2 (m 2 - m 3 ) 



ml log(m 2 ) + m\ log(r?7,2) {m\ + m\ — m\) 



log^ 2 ) m\(m\ — m\ — m\) + m\{—m\ + m\ — wij) | , (30) 



R 2 (ml, mj, ml) 



d 



dm\ 



[m\, m\, 777,3) = ( m 3 — m i ~~ m\)V^ (ml, m\, m 



+ (7773 — m 2 x — Trig) ^(m^m^m 



1 



+ —\ —ml(ml — m\) — m\{m\ — m\) — 2m\m\ log(m^) log(m2) 



+ 



m\ \og(m\) (l - - log(m 2 )^ + m\ log(m^) ^1 - - log(m^)^ 



f 2 1 2 2\ 

[m 1 + m 2 — 777 3 J 



+ log(m 3 ) ^1 — - log(m 3 )J m l (m l — m 2 — m 3 ) + m 2 (—m l + m 2 — m 3 ) 
— log^g) m\{m\ — m\ — 7773) log(m^) + m\(—m\ + m\ — m\) log(77?| 



Similarly, in the notation of Eq. (^) the coefficients of the expansion of Eq. fl27f) are 

1 



(31) 



2 \m\,m\,{nii + m 2 ) 2 ) 



l \m{, ml, (mi + 7712) 2 ) 
1 

~ 4 

V (0) (777 2 ,7772, (m X + m 2 ) 2 ) 

3 



-(m\m 2 + m x + m 2 ) 

-(777 1 777 2 + 777^ + 7T7 2 ,) 



m\ log(mi) + m l log( m 2) + ( m i + m 2) 2 log(mi + 777.2) 



7_ 
16 



(m\m 2 + m\ + m 2 ) 



+ 8 

1 



777 1 log(777i) + 777 2 log(7772) + (t77i + 7772) log(777i + 7712) 



^ [mi log(mi) - m 2 log(m 2 )] - -(m x + m 2 ) 2 \og 2 (m 1 + 777 2 ) 
o o 

1 



4 



[mi log(mi) + 777 2 log(m 2 )] (m x + 777 2 ) log(mi + 777 2 ) . 



(32) 



When solved recursively, starting from j = —2, the equations ( 29[ , |3~0 ,|3~I) are all of the 
form 



R 2 (ml, ml, ml) -^V^ (ml, ml, ml) 



(ml — ml — ml)V^\ml, ml, ml) 

(ml, ml, m 2 3 ) , (33) 
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where the g^\m\, m\, mj), j = —2, —1, 0, to be identified by comparison with the original 
equations ( p9| , |30y3TD , can be considered as known. 

Eq.(|33|) can in principle be solved by the quadrature formula 

V^\m\, ml, z) 



V^\m\,ml,ml) = R(m\,m\,m 



3; \ R(mlmlz) 

+ n mr ^ 2 , g {j) irnl,ml,x) \ , (34) 
Jz R A [m{,m 2 ,x) J 

where z is any initial point. Our initial point is (mi + m 2 ) 2 ; as for that value of z the 
two terms in the r.h.s. diverge, it is convenient to modify the formula by carrying out 
an integration by parts and then to take the z = (mi + m 2 ) 2 limit. When that is done 
Eq . (|34]) becomes 

if m ^ 2 2 n 

V^\m x ,m 2 ,m z ) =^,t [ R(m 1 ,m 2 ,m 3 ) J dx ^K,m 2) x) 

(m.i+m 2 ) 2 

+ (m 2 + m 2 -ml)g {] \m\,ml,ml)^ ; (35) 

that last formula, as can be easily verified, satisfies the differential equations ( f29| , |30l , |31~|) , 
with the initial conditions Eq. (|32|) , implied by the regularity at m.3 = m\ + m 2 . 

The integration can be actually carried out by using the change of variable 

2 o t 2 + 1 m\m 2 ( mA / m 2 \ 
x = m 1 + m 2 + m^mi = 1 1 H I I t H I , (36) 



t \ m 2 J \ m 1 



which can be inverted into 

t = t(x; ml, ml) 



x — (mi — m 2 y + \ x — (mi + m 2 



x — (mi — m 2 ) 2 — Jx — (mi + m 2 ) 2 

= (x — m? — ml + R(x, m?, ml)] . (37) 

2m\m 2 v ' 

The solutions of the equations (^,^,^l|) with the initial conditions given by the equations 
([32] ) are then respectively 

V { - 2 \m\,ml,ml) = --(m 2 + m 2 , + m 2 ) , (38) 
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\r(— 11/ 2 2 

l/ v ' [m 1 ,m 2 ,m 



1 f 3 . o 9 2\ 



m 2 log(mf) + wij log^ 2 .) + rn.g log(m 3 



(39) 



V (0 \ 



2 2 2^ 

m 1 ,m 2 ,m 3 ) 



1 
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Til 2 2 2\ r / 2 2 2\ 

R{m 1 ,m 2 , m i )L[m 1 ,m 2 , m 3 ) 



/ 2 i 2 , 2\ 

+ m 2 + m 3i ) 



+3 



m 2 log(m 2 ) + m\ log(m 2 ) + m 3 log(m 2 



+ 



m\ log 2 (m 2 ) + m 2 log 2 (m 2 ) + m 2 log 2 (m 3 ) 

(— m\ — m\ + m 3 ) log(m 2 ) log(m 2 ) 
+ (— m\ + m 2 — m 2 ) log(m 2 ) log(m 2 ) 
+ (m 2 - m 2 - ml) log(m 2 ) log(m 2 ) 



(40) 



where 



L(ml, ml, ml) = U 2 {-t 3 m 2 /m 1 ) + Li 2 (-t 3 m 1 /m 2 ) + ((2) + - log 2 (t 3 ) 



1 

+ 2 



log(t 3 + m 2 /mi) - log(t 3 + mi/m 2 ) + - log(m 2 /m 2 ) 



log(m 2 /m 2 ) , (41) 



with t 3 = t{m\] m\ , m 2 ), see Eq.(p7|) . It is to be noted that Eq.s (p8| - ^0D are expected 
to be symmetrical in the exchange of m 2 , m\, m 2 ; the symmetry, which is manifest for 
most terms, holds also for L(mf, m\, mj), Eq.(^TJ), even if in this case the symmetry is not 
manifest. Note also that Eq. (|4T|) remains valid in its present form for any real positive 
value of mi, m 2 and m 3 , provided that t 3 is given by the second of Eq.s (|37|) . In the 
limiting case mi,m 2 — > 0, finally, independently of the order of the limits one finds 



V (0) (0,0, 



m 



m 
1~6 



2 r 



7 

C(2) + - - 3 log(m 2 ) + log 2 (m 2 



(42) 



Analytical values of V(n, ml, m\, ml) were already presented in H, where a differential 
equations method based on mass derivatives was used, and in , where the results were 
obtained by using a Mellin-Barnes representation and its subsequent expansion in n — ► 4 
limit. Our results agree numerically with those presented in || in terms of Lobachevskiy's 
function and a related one; to check that it agrees also with the analytic formula of Jf) it 
is necessary to use the following relation between dilogarithms of different arguments 



Li 2 (-t 2 m 1 /m 3 ) + \A 2 {-t 1 m 2 / m 3 ) - log(-t 2 mi/m 3 ) log(-tim 2 /m 3 ) 
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I I 

U 2 (-t 3 mi/m 2 ) + U 2 (-t 3 m 2 /m 1 ) + 2((2) + - log 2 (t 3 ) - - \og{m\/ml) \og(m 2 2 /m 2 3 ) 

1 3 
+ -[\og(t 3 + m 2 /m 1 ) -\og(t 3 + mi/ m 2 )} \og(my mj) + -log 2 (mj/ mj) , (43) 

Z o 

valid for m 3 > mi +m 2 , with t± = t(m\; m\, m§) and t 2 = t(m^]mf,m^), whose validity 
can be established by standard methods. 



4 The expansion at p 2 ~ 0. 

As p 2 = is a regular point for the 4 functions Fi(n,m 2 ,m 2 ,m 3 ,p 2 ), i = 0, 1,2,3, they 
can be expanded in powers of p 2 around that point 

oo 

Fi^mlmlmlp 2 ) = ]T F i>k (n, ml, m\, m 2 3 )(p 2 ) k . (44) 

fc=0 



The very first term, according to Eq. (f24"|) is 

F , (n, ml, m\, m\) = F (n, m\, m\, m\, 0) = V(n, m 2 x , m\, m 2 ) , (45) 

where the vacuum amplitude V{n, m\, m\, ml) discussed and evaluated in the previous 
Section, Eq.fl2|) and Eq.s(|||g(]). Once V{n,m 2 ,m 2 ,m 2 ) is known, Eq.(g§ and Eq.@ 
can be combined to give, if {i, j, k) is any permutation of the indices (1, 2, 3) 



F ifi {n,ml,m 2 2 ,m 2 3 ) = Fi{n,ml,m 2 2 ,m 2 3 ,Q) 

1 



(n — 3)(m 2 — m 2 — m 2 k )V(n, ml, m\, m 2 ) 



R 2 (m 2 , m\, m 3 ) 

mj -m 2 + m 2 

+ \ n ~ 2 ) 7T2 T{n,m,)T{n, mj 

2mf J 

m 2 + m 2 -m 2 k 2 2 

+ ~ ' 2m^ T ( n ' m fc) T( - n ' m i ) 



— (n — 2)T(n, m 2 )T(n, m 2 k ) 



(46) 



showing that the three quantities F ii0 (n, mf , m\, ml), % = 1,2,3, are all expressed alge- 
braically in terms of V(n, ml, m\, ml) only. 

When the expansion Eq.(|4]) is inserted in the Eq.s (§) and (^), by imposing the equal- 
ity of the coefficients of the various powers of p 2 one obtains a set of equations for the quan- 
tities Fi^[n,ml,m\,m\). For k = one recovers of course the equations discussed and 
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solved in the previous Section; for k > one gets a set of linear algebraic equations for the 
F it k(n,m±,m2,rnl), which can be solved recursively in terms of the Fi t k-i{ n , m 1, m \i m V)- 
Therefore everything can be expressed in terms of the Fifi(n,m\,m\,rn%), which in turn 
depend, as already remarked, on the single vacuum amplitude V r (n, mf, m^mf); as a 
consequence, all the ^^(n, mf , m|, m|) can be expressed in terms of V(n, m\, m 2 ., m|) 
only 

The first iteration of the just sketched procedure gives 

n — 3 



F ,i(n,m v m 2 ,m 3 ) 



nR 4 (rrii, mf, m|) 



(n — 4) ( m 1 (m 2 + m 3 — mf) 



+ m\{m\ — m\ + 777.3) + m 3( m i + m 2 — m 3)) — 2(n — %}m\m\m\ 



V{n, ml,ml,ml) 



+ 



+ 



+ 



n-2 



nR 4 (rrii, m\, m|) 

n-2 
nR 4 (m 2 , ml, mf) 

n-2 
nR 4 (m 2 , m|, m§) 



(77 — 4) (7/73 + — 7772) (7773 — + m\) 

+2m\{m\ + wij — m\) T(n,m1)T(n,ml) 

[n — 4){m\ + m\ — m\){m\ — m\ + 777,3) 

+27772(777^ — 7772 + 7773) T(n, m\)T{n, m\) 

(77 — 4) (777^ + 7772 — Trig) (777^ — 7772 + 777,3) 

+2m\{m\ — m\ + m^) T(n, 7/72)7(77, 7/73) . (47) 



As already done for V{n,m\,m\,m^) ^ n Eq.(P%D, the F it k(n,m\,m\,m\) can be in turn 
expanded in (n — 4) as 



77 / 2 2 2\ 



C 2 {n)< 



Fk 2) ( m l,ml,ml) + 



+F$(ml,ml,ml) + 0(n-4) 



(48) 



where C(77) is given as usual by Eq.(0) and Eq. ([□]). From Eq. (|47|) and Eq.s ([38|440D one 
obtains 

F^ 2 \m 2 1 ,m 2 2 ,ml) = , (49) 
F^ 1) (m 2 1 ,m 2 2 ,ml) = ^ , (50) 



32 



z+o)/ 2 2 



1 f 4777?r77ar77o 

J - •' -L{ml, mlml) 



32 I R 3 (mi, m\ 



3 
4 
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R 2 (m 2 , 771%, ml 



2/2 2 2\ l / 2\ 

m 1 {m 1 — m 2 — m 3 ) log{m 1 ) 



+m\{—m\ + m 



log(m 



l)+rn\{ 



-ml — ml + ml) log(m 



(51) 



This result is in numerical agreement with the first term in the p 2 /ml expansion in M. 
It is straightforward to obtain, when required, any other term of the double expansion in 
p 2 and (n — 4) . 



5 The expansion in (n — 4) for arbitrary p 2 . 

Quite in general, by using the same argument leading to the expansion of Eq.(^), the 
functions Fi(n,m 2 ,ml,ml,p 2 ),i = 0, 1,2,3, are expected to develop up to a double pole 
in (n — 4) and can be expanded, also for arbitrary p 2 , in the by now customary form 

IP/ 2 2 2 2\ ri2i \ \ 1 Z7(- 2 )/ 2 2 2 2\ 

F i (n,m 1 ,m 2l m 3 ,p ) = C (n)| ^ '{m 1 ,m 2 ,m 3 ,p ) 

+ j^j^W^^Ip 2 ) + if ViX^P 2 ) + 0(„_4)j (52) 

When all the terms of Eq.s(||, 0) are consistently expanded in (n — 4), one obtains a set 
of equations for the residua and the finite part in (n — 4). Those equations can then 
be solved recursively, starting from the double pole. The equations corresponding to the 
poles in (n — 4) will not be written explicitly here for the sake of brevity; it is however 
easy to verify that the solutions of those equations, corresponding to the coefficients of 
the poles in (n — 4) of Fo(n, m 2 , m 2 ,, m|,p 2 ) are 

rp(-2) / 2 2 2 2\ ^ / 2 , 2, 2\ 

F V J (m 1 ,m 2 ,m 3 ,p ) = - + m 2 + m 3 ) , 

■j—. ( — 1) / 9 2 2 2\ I ^ ^ / 2 2 2\ 

^0 (™!,™2> m 3>:P ) = gS J + +m 2 + m 3) 

m 2 log(m 2 ) + log(m 2 ) + ml log(mg) j , (53) 

from which the corresponding coefficients (mf , m|, ?n|,p 2 ), for i = 1,2,3 and = 
—2, —1 are immediately obtained by using the Eq.s(|3]), 

Ei(-2)/ 2 2 2 2\ ^ 

^ (-1) (m?,TO2,mi,p 2 ) = -^ + ilog(m 2 ). (54) 
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The above functions, indeed, satisfy the differential equations and the initial condition at 
p 2 = provided by Eq.s( p^ - |^0| ) , therefore they are the exact solutions for arbitrary p 2 . 

On account of the above results for the singular parts, from Eq.s([|, [1]) one obtains for 
the finite parts F^ (m 2 , m 2 ,, m\,p 2 ) of the 4 amplitudes the equations 

, d 



t7t(0)/ 2 2 2 2 
P Q^ F ( m D m 2) m 3>P 

,2 77(0), 



itt(O)/ 2 2 2 2 \ i 77t(-1)/ 2 2 2 2\ 

FJ '{m 1 ,m 2 ,m 3 ,p ) + F V ; (m 1 ,m 2 ,m 3 ,p ) 



' ; 2 2 2 2., 277(0)/ 2 2 2 2, { ;>!"), 2 2 2 2, / r- r , 

+ 777^ '(m 1 ,m 2 ,m 3 ,p ) + m 2 F 2 v ; (m 1 ,m 2 ,m 3 ,p ) + m 3 F 3 v ' [m x ,m 2 ,m 3 ,p ) , (55) 



277(0), 



and 



.j > > .) j (9 ifO'l / *)\ 

8D(m 1: m 2 ,m 3 ,p )p -^F> (m^m^m^p ) 



at\! 2 2 2 2\77i(— 1)/ 2 2 2 2 



1 72> / 2 2 2 2 

+P i)0 (m 1 ,m 2 , 777,3, p 



lQF$ ) \m\ 1 mlmlp 2 ) + 28F ( 1} (m?, m^, m^,p 2 

+12F ( - 2) K,m2,m2,p 2 ; 



+8P M (m 2 ,m^,m^,jj 2 ) f/ ^ (m 2 , m 2 , m 3 , p 2 ) + Fj l ' [mi,m 2 ,m 3 ,p~ 



+8Fj(m 2 ,m 2 ,m 3 ,jj 2 ) 
+8F jfc (m 2 ,m 2 ,m 3 ,77 2 ) 



'77(0)/ 2 2 2 2 \ 1 77(-l)/ 2 2 2 2 

p- (m^ m 2 , m 3 ,p ) + t > (m-L, m 2 , m 3 , 77 



77(0)/ 2 2 2 2 \ 1 uif-l)/ 2 2 2 2\ 
2_2 1. ,2 



i 1} (mj, 7775,7715, jr 



+Qi, i (m 1 ,m 2 ,m 3 ,p ) ra^m*. [log(m i ) + log(?77 fc 
+Qij(ml,ml,m 2 3 ,p 2 ) m 2 m\ [log(m 2 ) + log(m^ 
+Q i ,k(m 2 1 ,ml,ml,p 2 ) m 2 m 2 log(m 2 ) + log(m 2 



(56) 



where fc) is any permutation of (1,2,3) and P^m 2 , m|, m^p 2 ), Q i j(rn'l,rn 2 ,'ml,p 2 
are the polynomials of Eq. (^) . 

In the notation of Eq.(|52|), for j = — 



1,0 



Fo^^m 2 , 7772, 7773,0) = \/ (i) (m 2 , 777 2 , 777 



2 1 
3 J J 



(57) 



where the V^(mf, m 2 , m 3 ) are given in Eq.s( |3"8"| - f4"0| ), so that FQ°\ml, m 2 , m 2 , 0) is known; 
similarly, the Ff\m\, m 2 ,, m 2 , 0), z = 1,2,3, can be obtained by using either Eq . (|25|) or 
Eq.([46|) (properly expanded in the 77 — > 4 limit). As a conclusion, Eq.(|55|) and Eq.(|56|), 
together with the initial conditions at p 2 = just discussed, are free by construction from 
the singularities of the n — > 4 limit and therefore ready to be used for the numerical 
evaluation of the finite parts of the sunrise master integrals. 
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For convenience of later use, we introduce one more quantity jF (m 2 , mf, m 2 ,p 2 ), de- 
fined by 



J r o{m 1 ,m 2 ,m 3 ,p ) = hm F (n,m 1 ,m 2 ,m 3 ,p 

n— >4 L 



2 2 2 2\ 

2 _2 2\ 2 77 /„ _2 2 2n 



- F 0)0 (n, m 1 ,m 2 , m 3 ) - p F ,i(n, m l5 m 2) m 3 )J , (58) 

where F 0i o(w, mj, m 2 , mj) and Fo,i(^, mf, m|, mf) are given by Eq. (ff5j) and Eq. P^) . 

It is to be noted that the terms of Fo(n, mf , m|, m|, p 2 ) singular in (n — 4) are at most 
linear in p 2 , as shown in Eq.s(|53"D, and therefore entirely encapsulated in the combination 
F 0t0 (n, m 2 , m 2 ,, m|) + p 2 F 0>1 (n,m 2 ,m 2 ,ml); for that reason J r (m'l,'m 2 ,rnl,p 2 ) is finite in 
the n — > 4 limit. 

From the very definition Eq. fl58|) , jF (m 2 , m|, m|,p 2 ) has a double zero at p 2 = 0; 
therefore it satisfies the doubly subtracted dispersion relation 

T^m\,m\,m\,p 2 ) = (p 2 ) 2 f — - ImF (4,m^m^m^-«) , (59) 

(mi+m2+m3) 2 

which is obviously finite at n = 4. 



6 The analytic calculation for m 2 = m 3 = 0. 

We discuss now the case in which all the masses but one are equal to zero; the remark- 
ably simpler equations obtained in that limit can be solved analytically, providing useful 
information for the study of the large p 2 behaviour, even for arbitrary masses, which will 
be worked out in the next Section. 

When putting m 2 = = 0, and mi = m for simplicity, the system of the 4 equations 
Eq.(||) and Eq.(^) collapses into just two equations: 

p 2 ^jF (n, m 2 , 0, 0, p 2 ) = m 2 F 1 (n, m 2 , 0, 0, p 2 ) + (n - 3)F (n, m 2 , 0, 0, p 2 ) 
op z 

p 2 (p 2 + m 2 ) ■^F 1 (n, m 2 , 0, 0, p 2 ) = -\ (n - 3) (3n - 8) F (n, m 2 , 0, 0, p 2 ) 



+ 



2p 2 (n - 3) - (p 2 + m 2 ) (-n - A^j 



Fi(n,m 2 ,0,0, p 2 ) . (60) 



The system so obtained is homogeneous, and as such it cannot provide information on 
the initial values; but the required initial conditions are still supplied by the proper 
n^2 = n^3 = limit of the results of the previous sections. 



17 



To solve the system, we observe that it is of course equivalent to a second order 
differential equation in a single function. If F (n, m 2 , 0, 0,p 2 ) is kept as the independent 
function, eliminating F 1 (n, m 2 , 0, 0,p 2 ) from Eq.s(|60|) one obtains 

d 2 f ti \ 

p 2 {p 2 + m 2 )——F {n, m 2 , 0, 0,p 2 ) + -2p 2 {n - 3) + -{p 2 + m 2 ) ^-zF {n, m 2 , 0, 0,p 2 ) 



d(p 2 ) 2 
(n-4)(n-3) 



dp 2 



F (n,m 2 ,0,0,p 2 ) = . 



(61) 



By introducing jF (m 2 , 0, 0,p 2 ) according to Eq. (|58D , with the values of F 0t0 (n, 0, 0, m 2 



and Fo^(n, 0, 0, m 2 ) corresponding to Eq.s(flT)|),(f47|) in the mi = vri2 = limit given by 
Eq. (p8|) and Eq. ([22|) , one has 



F (n,m 2 ,0,0,p 2 ) 



1 

+ 16 



P 



C\n) 

3 !i 

-- + -log(m 



m 



+ 7 — ~ — TTT^T (p 2 + 6m 2 - 4m 2 log! 



(n-A) 2 8 (n-4)32 
:! ^ - m 2 (((2) + - - 3 log(m 2 ) + log 2 (m 2 ) 



m 



+^o(m 2 ,0,0,p 2 ) + 



(62) 



In terms of jF (m 2 , 0, 0,p 2 ) Eq. (|6T|) then reads 

?V + m2 )^ji^o( m2 > °> °> P 2 ) + 2m2 ^^o(m 2 , 0, 0, p 2 ) = ^p 2 , (63) 

which is in fact a first order differential equation in -J=^ Tq{vp? ^ 0, 0,p 2 ). The values of 
^(m 2 , 0, 0,p 2 ) and its derivative at p 2 = are by construction 



Jo(m 2 , 0,0,0) = 0, 
^o(™ 2 , 0,0,0) = 0. 



(64) 



With the substitution 



d_ 
dp 2 ' 



— Jo(m 2 ,0,0,p 2 ) 



(p 2 + m 2 ) 2 



-G(m 2 ,p 2 ' 



( P 2 y 



(65) 



Eq.(^) becomes 



_d_ 

dp 2 



G(m 2 ,p 2 ) 



1 (P 2 ) 2 
32 (p 2 + m 2 ) 



(66) 
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The integration of Eq . (|66"D is trivial; the result in terms of jF (m 2 , 0, 0,p 2 ), accounting for 
the second of the initial condition of Eq.(^), reads 



!^ K,o,oy) = J_ 



(p 2 + m 2 ) 2 , p 2 . m 2 3 

7 2X2 M 1 + —)- — - o 

[p ) vn A p z 2 



(67) 



One more integration and the first initial condition of Eq.(^) give 
Mm 2 A0,p 2 ) = lJ-2m 2 U 2 (-Q + (p 2 - log (l + Q - 5 -p 2 + m 2 \ . (68) 



32 



The same result can be obtained, by means of the doubly subtracted dispersion rela- 
tions Eq.(|59|) for the specific values of the masses, 



Jo(m 2 ,0,0,p 2 ) = (p 2 ) 2 ( dU - ImF (4,m 2 ,0,0,-n 



U A {U+p Z J IT 



(69) 



The Cutkosky-Veltman rule gives for the imaginary part 



1 , 9 , 1 r „ u — b b — m 2 

- ImFo 4,m 2 ,0,0,-w = — db — 

7i 16 J 2 u b 



1 I 2 

1 \ u — m 



u 



m !og — ? ? ; 



16 [ 2m 

substituting in Eq. (|69|) and carrying out a last integration Eq.(|68|) is recovered. 



(70) 



7 Large p 2 expansion. 

By using the differential equations (^|J^) it is also possible to study the asymptotic expan- 
sion of the 4 functions Fi(n, m 2 , m|, m 2 ,p 2 ) for large values of p 2 . As p 2 = oo is a singular 
point of the equations, we look for a tentative asymptotic expansion of Fq in the form 

OO -| 

Ft\n,mlmlmlp 2 ) = (p 2 ) a £ ^(n, m 2 ,m 2 ,m 2 )-i- . (71) 

Due to the linearity of the equations, the actual solution is built by summing a suitable 
linear combination of all the possible expansions. The corresponding expansions of the 
other 3 functions F^ a \n } m|, m 2 , m 2 ,p 2 ), i = 1,2,3, can be obtained through Eq.(|). 
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By inserting the expansion Eq. (|7lD into Eq.@ and Eq.(|7|)), one obtains a set of 4 
equations, which must be satisfied, for each separate power of p 2 , by the coefficients of 
the expansion. 

A first solution corresponds to a = —1; we write it as 



It corresponds to the socalled regular solution, hence the superscript (r) in <p[ (we write 

(r) (r) 

(f) k rather than § k for uniformity with later use). Its coefficients are completely de- 
termined by the non-homogeneous part of the equations; correspondingly, all the terms 
behave for large p 2 as strictly integer powers of 1/p 2 , which is exactly the behaviour of 
the various terms occurring in the asymptotic expansion of the non-homogeneous part. 
The first coefficient of the expansion is 

(p^\n, m\, m\, ml) = T(n, m 2 )T(n, m^) + T{n, ml)T(n, ml) + T(n, m 2 )T(n, ml) ; (73) 

(r) 

all the other coefficients <j) k can be similarly obtained, when needed, by following the 
procedure just outlined. 

Other solutions correspond to different values of a; if a is not equal to —1, the leading 
coefficient must satisfy the homogeneous system of equations 




(72) 



H m ?^~^ $ o {n, ml, ml, ml) + (a - n + 3) % a '(n, ml, m\, m 2 3 ) = 




(74) 



There are two solutions of the system ([74]), corresponding to the values a — (n 
a = (n — 4)/2. Correspondingly, Eq.s (|74]) become 



3) and 




(75) 



and 




(76) 
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the above values of a are also called "fractionary", and the corresponding solutions sin- 
gular, hence the superscripts (si), (s 2 ), as opposed to the superscript (r) for the "regular" 
solution corresponding to the integer value a = — 1 of Eq. ([72]) . Due to the homogeneity 
of the equations for the coefficients of the "fractional" powers, the two leading coeffi- 
cients $Q l \n, mf, m 2 , m|) cannot be determined by the equations ( [75]) and Q75D alone; 
nevertheless, the equations can still provide much valuable information. 

Eq.(|75f) shows indeed that $g 1 \ n i m l m l m i) * s independent of the masses 

$ { Sl) (n,mlml } ml) = $ (si) (n) , (77) 

so that the still unknown function & si \n) depends only on n. 

The case of Eq. ( [76|) is slightly more complicated; by investigating the condition on the 
leading (p 2 )( n ~ 4 )/ 2 term implied by Eq.fl9~5|) of the Appendix, whose l.h.s. is the second 
derivative of F (n, m 2 , m 2 , m|,p 2 ) with respect to m 2 and mf, one finds the relations 

$( S2) (n,m 2 ,m 2 ,m 2 ) = (i, j = 1, 2, 3; z# j) . (78) 



dmfdm'j 



According to Eq.([75]), the solution of Eq.([7j]) can be written as 

<!> { S2 \n,ml,m 2 2 ,ml) = $M(n)$ 2 \n, m 2 , m 2 , m 2 ) , 

^(n, ml ml ml) = m\~ 2 + m n 2 ~ 2 + m n 3 ~ 2 , (79) 

where & S2 \n) is a function of n only. 

For each "fractional" power solution, all the next-to-leading terms of the expansion 
are determined in terms of the leading terms, by solving recursively the systems of linear 
equations; it is therefore convenient to factorize the (as yet arbitrary) function & Si >(n), 
defining 

$jf } (n, ml ml ml) = ^ Si \n)^ i \n,mlmlml) . (80) 
An explicit calculation then gives for the first few terms (which we will use later) 

i(st)/ 2 2 2\ \ n ~ 3)(3n — 8) / o 2 2\ 
01 {n^m^m^m^ = ( n _ 4 ) [rn 1 + m 2 + m 3 ) , 

TDl 2 2 2\ i d2/ 2 2 2\ 

B(m 1 ,m 2 ,m 3 ) + R (m 1 ,m 2 ,m 3 ) 



.(-0/ 2 2 2n (w-3)(3n- 
2 '(n,m 1 ,m 2 ,m 3 ) = _ 



(n-4)" 



^KXX) + §(„ -4) + ^(m 2 ,m 2 ,m 2 )) 



(81) 
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where 
and 



B (ml 



,m\,m\ 



A ( 2 2 i 2 2 i 

4 I m l m 2 + m l m 3 + 



m\m\ } . 



(82) 



J.{ s 2) f 2 2 2, 

01 '(n,m 1 ,m 2) m 3 ) 



n 



-(3n — 8) ( m l + m 2 + m 3 J <% (n,m 1 ,m 2 ,m 



+ 2(n - 2) (m^ 2 (m\ + 2/7*2 + 2mg) + m^ 2 (2m] + + 2mg) 



+ ml-\2m{ + 2ml + ml) 



(83) 



Any subsequent term can be easily obtained, when needed, by the same procedure. 

As for arbitrary n the three solutions satisfy the equations independently of each other, 
the actual asymptotic expansion of Fo(n, m 2 , m|, m 2 ,p 2 ) can be written as their sum 

1 \ 



t-.(oo)/ 2 2 2 2^ 

F$ J {n,m 1 ,m 2 ,m 3 ,p ) 



P 



(/r 4 $ (si) w 1 + E 



k=l 



ds2). 



n,m\,m\,ml)j^ 

,2 _2 „2\ 1 



fc=0 / 
1 00 1 



fc=0 



(p 



2U 



J4) 



the above expression depends on the two functions of n, $( si )(n) and $( S2 )(n), which are as 
yet undetermined. To fix them, we must supply more information on F Q (n, m 2 , m 2 , m 2 , p 2 ). 
As the Fi(n,m 2 ,m 2 ,ml,p 2 ),i = 0,1,2,3, were studied in the n — > 4 limit, Eq.(|2|), we 
expand also Eq.flS4l) around n = 4, writing in particular & Sl \n) and & S2 \n) as 



$<-*>( n ) = C 2 {n) ]T(n- 4) fc $ ( ^ } (i = 1, 2) ; 



J5) 



l=-2 



as always C(n) is not expanded, appearing as an overall factor. 

By comparison with the coefficients of the singularities in (n— 4), which are exactly known 
and given in Eq.([53|), the following values are found 



(J>0l,-2) = q ^ $(s 2 -2) = _. 



$0i -1) 



32 



$0»2,-i) _ _4$(si,o) 



32 



i6) 



Unfortunately some additional information has still to be supplied, besides that coming 
from the singularities in (n — 4), for fixing the other as yet unknown constants $( Sl '°), 
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$(s2,o) anc j $(si,i). the i as ^ quantity is the coefficient of the term in (n — 4) in Eg. (|8~5"D, 
which is also needed due to the presence of the 1 / (n — 4) factor in Eq. (j5TD . 

As those constants are independent of the masses, they can be determined through the 
comparison with the asymptotic expansion of a simpler case, namely F (n, m 2 , 0, 0,p 2 ), 
discussed in the previous Section, whose explicit analytic value is given in Eq.flBSp. Its 
large p 2 expansion reads 

Ft\n,m 2 Mp 2 ) = C^{n)[-- r ~^ + T^TT^ ft + ^ " ^ bg(m2 



1 

+ 16 



(n-A) 2 8 (n-4)32 

P 2 ( + \ lo g(P 2 )) + 1o S 2 (p 2 ) - ™ 2 log(m 2 ) logQo 2 ) 

1 5 
+ m 2 ( — - log 2 (m 2 ) + 31og(m 2 ) - - 



- ^2 ^(P 2 ) _ !og(> 2 ) + ~ 



+ ... . (87) 



Let us recall that, when using Eq. (|84|) and Eq. fl35|) , the logarithmic factors are generated 
in the n — > 4 limit through expansions of the kind 

' rCP 2 )*"" 4 = 7-^7T? (l + (»» - 4 ) ^(P 2 ) + i(n - 4) 2 logV) + • 



(n-4) 2V ^ y (n-4) 2 V 2 

By comparing Eq.(]87f) and the n — > 4 limit of Eq.(jS4|), when two masses are taken equal 



to zero, Eq. (|86D is completed into 



$(»l,-2) = ^ $(«2 -2) = _i 

$(•1.-1) = _L ; $(*2 -i) = I 
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$(•1,0) = _ J£ $ ( S2 ,o) , 4 $( S i,i) = / 88 a 

128 " 128 ' y ! 

It is apparent from the above equations that only a particular combination of $( S2 >°) 
and $( Sl>1 ) is determined by the comparison, but only that combination is needed for 
expressing Fq(ti, m\ , m 2 ,, m^,p 2 ) up to the constant term in (n — 4) ( i.e. the separate 
knowledge of the two terms is not needed up to the degree of expansion in (n — 4) and 
p 2 , considered here). 

Summing up, for arbitrary masses and in the n — > 4 limit, up to terms 1/p 2 included, 
the asymptotic expansion Eq. fl84]) can be written as 



Ft\mlmlmlp 2 ) = C\n)L^Ft^ 2) {ml m 2 , m 2 ,p 2 ) 
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1 1 
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(89) 

(90) 
(91) 



(92) 



-3B(ml ml m 2 ) - l(LR 2 (m 2 , m 2 , m 2 ) + 4(m 2 - 2m 2 - 2m 2 )m 2 log(m 2 ; 



+ 4(m 2 — 2m 2 — 2m 2 )m 2 log(m 2 ) + 4(m 2 — 2m 2 — 2m 2 )m 2 log(m 2 
+ 8m 2 m 2 log(m 2 ) log(m 2 ) + 8m 2 m 2 log(m 2 ) log(m 2 ) 
+ 8m 2 m 2 log(m 2 ) log(m 2 



(93) 



The coefficients of l/{n — 4) 2 and l/(n — 4), Eq. (^0) and Eq.(^TJ), are of course identical to 
the exact ones in Eq. fl53f) ), as their expansions in 1/p 2 contain only a linear and a constant 
term. The results given in Eq.s (|92|j93|) agree with the first terms of the corresponding 
expansion given in |§. 

Higher order terms, when needed, can be added with little effort by following the 
above discussed procedure. 
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Appendix. 

We list here a few ready-to-use formulae, obtained by solving the integration by part 
identities for the amplitudes of the sunrise self- mass graph (see also ||). The notation 
for the amplitudes is given in Eq. ([!]). Further relations can be found by permutation of 
masses. 



n/2 2 2 2\ A I 2 2 22 

D{m x ,m 2 ,m 3 ,p )A{n,m x ,m 2 ,m 3 ,p 

-D(ml, ml, m 3 ,p 2 )Fi(n, ml, m\, m\,p 2 ) 



-3,-1, 



-1,0,0) 



4m 2 



+ (n — 3)(p 2 + m\) 4(m 1 m 2 + m\m\ — m 2 m 3 ) 

/ 2 i 2 , 2 , 2\2l 77 I 2 2 2 2\ 

— [p + m x + m 2 + m 3 ) t 1 {n,m 1 ,m 2 ,m 2) ,p ) 



mr, 



+ (n — 3)(p + m 2 ) — 2 4(m 1 m 2 — m 1 m 3 + m 2 m 3 ) 



/ 2 i 2 i 2 i 2\2 rrW 2 2 2 2\ 

[p + m l + m 2 + m 3 ) r 2 {n,m l ,m 2 ,m 3 ,p ) 



( o \ / 2 2\ T^/ 22 22 22'\ 

+ [n — 6)(p + m 3 ) — 2 4(— m x m 2 + m 1 m 3 + m 2 m 3 ) 



in 



+ (n-3)(3n-8)-l r 



/ 2 i 2 i 2 i 2\2l rrW 2 2 2 2\ 

— (p + m x + m 2 + m 3/ ) r 3 {n,m l ,m 2 ,m 3 ,p ) 

^ 2 2 2 i / 2 2 i 22, 2 2\/ 2 , 2 , 2 , 2\ 

-4m 1 m 2 m 3 + (m 1 m 2 + + m 2 m 3 ){p + m l + m 2 + m 3 ) 

1 / 2 , 2 , 2 , 2\3l 17 / 2 2 2 2\ 

- -(p +T7T,! +m 2 + m 3 ) \F (n,m 1 ,m 2 ,m 3 ,p ) 



+ 



+ 



+ 



{n-2f 



2 



2\2 



12 2\ I 2 2\ / 2 , 2, 2 2\ 

(m 3 - - m 2 ) - -(p + m 1 + m 2 - m 3 ) 



mf 
{n-2f 

2 

mf 



2 r 



<• 2 2\/ 2 2\ / 2 i 2 2 , 2\2 

v m 2 - - m 3 ) - -(p + m 1 - m 2 + m 3 ) 



T{n, m\)T{n, m 2 
T{n, m\)T{n, ml 



(m l — m 2 )(m 1 — m 3 ) (p — m 1 + m 2 + m 3 ) T(n,m 2 )T(n,m 3 ) , (94) 
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